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Two-particle tunneling in synchronous and asynchronous regimes is studied in the framework 
of dissipative quantum tunneling. We show that the use of the proposed model is justified by a 
comparison with realistic potential energy surfaces of porphyrin and experimental dependence of 
the reaction rate on temperature. The critical temperature Tc corresponding to a bifurcation of 
the underbarrier trajectory is determined. The effect of a heat bath local mode on the probability 
of two-dimensional tunneling transfer is also investigated. At certain values of the parameters, 
the degeneracy of antiparallel tunneling trajectories is important. Thus, four, six, twelve, etc., 
pairs of the trajectories should be taken into account (a cascade of bifurcations). For the parallel 
particle tunneling the bifurcation resembles phase transition of a first kind, while for the antiparallel 
transfer it behaves as second order phase transition. The proposed theory allows for the explanation 
of experimental data on quantum fluctuations in two-proton tunneling in porphyrins near the critical 
temperature. 

PACS numbers: 03.65 Xp, 03.65 Sq, 31.15 Gy, 31.15 Kb, 73.40 Gk, 82.20 Xr 



I. INTRODUCTION 



Quantum tunneling dynamics of a particle interacting with a heat bath is one of the important problems of modern 
condensed matter physics Q-jl^- The interest to this problem is related to studies in low-temperature superconductive 
tunnel junctions p, Ip, a dissipative quantum tunneling in crystals ^^'^ low-temperature chemical reactions 
P, "2?, '25', '2^, '2?, "28", '29', 'SCfl . In low dimensional systems an effective mass approximation often fails. Thus, quantum 
tunneling with dissipation becomes an important tool in the description of electron trans fer 13]. In many physical 
processes tunneling of two particles occurs. For example, Semenov and Dakhnovskii |2^ l29f found the bistability of 
the tunneling trajectories in a two-proton transfer. Later, Benderskii and co-workers extended their investigations to 
different two-dimensional potentials 30] . Some features of the two-dimensional tunneling dynamics were studied in 
interacting Josephson junctions [^. 

Another example of two-particle tunneling is a two-proton transfer in porphyrin systems l28l l29ll . Porphyrins 
are important molecules in Biology lU and a new area of Electronics - molecular wires and devices |32] . There are 
experimental data OH^H^H^] that clearly indicate that a tunneling instability occurs at some critical temperature. 
Such a behavior in the rate constant demonstrates the existence of bifurcations in two-dimensional underbarrier 



•Electronic address: aring azin@mail.kz| Also at Kazakhstan Division, Moscow State University, Moscow 119899, Russia. 
^URL: http://www.emu.kz 

Electronic address: yurid@uwy o.edu| 
§ Electronic address: physicsOdiamond. stup.ac.ru' 

^Electronic address: physics@diamond.stup.ac.ru Also at Institute for Basic Research, P.O. Box 1577, Palm Harbor, FL 34682, USA. 



2 



trajectories. Apparently, this effect requires a thorough investigation when the protons interact with a thermal bath. 

The effect of cleavage in the temperature dependence of the current for two-dimensional systems of interacting 
Josephson junctions was identified in Ref. 0. However, there was no experimental confirmation of this effect. In the 
present paper we identify and study such an effect of unstable character observed experimentally in tunneling constant 
as a function of temperature for reactions with porphyrin. Also, we predict a stable cleavage in the temperature 
dependence of current in the artificial two-dimensional system, a pair of weakly coupled charged quantum dots with 
single electrons tunneling to superlattice substrate, and the associated fluctuations of quantum origin. We study the 
emerging effect of chaotization near the critical point of unstable cleavage for porphyrin. 

In this paper, we continue to study bifurcation effects in two dimensional tunneling. In particular, we consider two- 
proton correlations of different types within the framework of dissipative quantum tunneling (instanton) approach. 
Moreover we discuss a fine structure of bifurcations for systems with various potential energy surfaces ,2^ ,2^ .2^ 

The layout of the paper is as follows. 

In Section II, wc introduce two-dimensional model potential energy surfaces for a pair of interacting particles. In 
Section III we present experimental aspects of the model. In Section IV, we study the effect of temperature on the 
tunneling rate. In Sections V and VI, we calculate the rate for the parallel and antiparallel tunneling and provide the 
analysis of the origin of bifurcation. The effect of a promoting mode (an environment) is studied in Section VII. 



II. TWO-DIMENSIONAL POTENTIAL ENERGY SURFACES 

Consider two charges that tunnel in two independent double- well potentials U{qi) and U{q2) presented as follows 

U{q,)^^LU^iq,+af9{-q,) 



-AI+^uj\q,-b)^ 



0{q^), * = 1,2, (1) 



where the sum a + b determines the length of a "link" in the corresponding macrocluster fragment; A/ — cj^(fo^ — a^)/2 
is a bias (an asymmetry parameter of the potential); 0{qi) is a step function, and lo is the frequency (see discussion 
m Ref. |27|). The mass is absorbed into the definition of o. 

The interaction between two charges, e.g., protons, is considered in a dipole-dipole approximation 

T4„t(gi,92) = -|(gi-g2)', (2) 

where a is a positive constant. We use the same interaction potential as in Ref. 0. 

Thus, the total two-dimensional potential energy surface for parallel tunneling normalized by w^, is given by 

Up{qi,q2) = — ^ = [qi + a) 0{-qi) 

+ [-{b" - a^) + (qi - bf] 0{qi) + {q^ + afOi-q^) 

+ [~{b' - a') + (92 - &)2] 9{q2) - ^(gi - q2)'. (3) 

Here, a* — 2a/uj'^ is the dimensionless parameter, a* < 1, a ~ e^/(ei?Q), e is the electron charge, Rq is the separation 
distance between the reaction coordinates qi and q2 of the tunneling particles; e is the dielectric constant. The form 
of potential energy surface Q is shown in Fig. ^ 

For antiparallel transfer, the two-dimensional potential energy with the interaction term can be defined as 

jj I X 2Ui,{qi,q2) 2af \ 

U^{qi,q2) = 5 = [qi + a) 0{~qi) 

+ [-{b^ - a^) + (gi - bf] e{q,) + {q2 - afe{q2) 
+ [-(62 - a2) + {q2 + bf] 9{-q2) -^(91" 92)^ (4) 

where a* = 2a/(M^ is a dimensionless parameter, a* < 1. The potential Q is depicted in Fig. [3 The main difference 
between the potential energy surfaces (O and Q is in the initial location of the second particle, ±a. 
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FIG. 1: Asymmetric potential energy surface Q for parallel tunneling; a — 2, b — 2.5, a* = 0.0001. A and B indicate the 
initial and final states of the particles, respectively. The minimum of the potential at B is lower than that of at A. The other 
two (intermediate) minima are lower than those of at A and higher than at B. 

The potential energy Q can be referred to as " parallel" while the potential energy Q as " antiparallel" . We also 
define a "symmetric" potential energy as a particular case of the potential Q under the condition a — b, i.e., 

Us{quq2) = 2 = Wi + ^^91 

+ {qi - afe{qi) + {q2 - af9{q2) + + afe{-q2) 

-\{qi-q2)\ (5) 

where the constant a** < 1. The potential energy © is presented in Fig. 13 

As shown in Refs. \^ , such model potential energy surfaces well describe the dynamics in porphyrins (see 

also Sec. Ill below). 

In many practical cases, the effect of a bath on particles' tunneling should be also included. In the next section we 
consider two particles embedded into a harmonic medium and linearly interacting with the bath modes. 

III. EXPERIMENTAL ASPECTS OF THE MODEL 

As it was mentioned in Introduction, adiabatic tunnel dissipative transfer appears to be very important in various 
physical and chemical processes [^-|30f. It is of much interest to consider two-particle tunnel dissipative transfer 
of coupled charges. For example, a two-proton transfer can be characterized by a simultaneous tunneling and the 
changing of electronic structure which forms and breaks chemical bonds 0; HE HOl- One can consider, e.g., 
transfer of two protons in the dimer of 7-azaindole jl^; see Fig. 01 

It is particularly interesting to consider two-dimensional tunnel dynamics of coupled protons in a symmetric case of 
porphyrin molecules. For this symmetric reaction an unstable interchanging of the regimes of tunneling (synchronous 
regime is replaced by asynchronous one as temperature decreases) was observed for the first time. The interchanging 
corresponds to the instability cleavage in the dependence of the rate constant on temperature. In some aspect this 
picture resembles the effect predicted for interacting Josephson junctions j^. In Fig. we schematically present the 
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FIG. 2: Asymmetric potential energy surface I^J for antiparallel tunneling; (a) a = 2, 6 = 2.3, a* — 0.1 (top panel); (b) a = 2, 
h — 2.3, a* = 0.5 (bottom panel). A and B indicate initial and final states of the particles. The minimum of the potential at 
B is lower than that of at A. The other two (intermediate) minima are higher than those of at A and B. 




FIG. 3: Symmetric potential energy surface @; (a) a — 2, b = 2, a** — 0.1 (top panel); (b) a = 2, b — 2, a** — 0.5 (bottom 
panel). A and B denote the initial and final states of the particles. The minimum of the potential energy at B is equal to that 
of at A. The other two (intermediate) minima are higher than those of at A and B. 
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FIG. 5: A view of two equilibrium configurations, I and II, for the porphyrin molecule with the transfer of a pair of hydrogen 
atoms between nitrogen atoms. The intermediate (nonequilibrium) configuration III is characterized by the hydrogen atoms 
residing on the neighbor nitrogen atoms. 



IDCO/T 



FIG. 6: Temperature dependence of the observed and predicted rate constants in the reaction of the internal transfer of 
hydrogen (or deuterium) atoms in porphyrin molecule. Open and filled circles denote experimental data for the hydrogen and 
deuterium atoms respectively 35]. The lines denoted by IH and ID represent predicted curves for a synchronous transfer 
of weakly coupled hydrogen and deuterium atoms respectively. The lines denoted 2H and 2D represent predicted curves for 
asynchronous transfer of weakly coupled hydrogen and deuterium atoms respectively. All the coupling constants J were scaled 
to reproduce the results for hydrogen atom transfer at T = 300 K. The used parameters are: the stretching and bending 
frequencies for the hydrogen atom are 3318 and 2482 cm~^, for the deuterium atoms are 1223 and 1094 cm~^ respectively l3^ : 
H and D stretching anharmonics and reduced masses are —80 and —50 cm~^, and 1 and 2 respectively; the skeletal frequency 
is 380 cm~^ and the reduced mass is 2.8. The equilibrium separation N— H (N— D) is 400 pm. Endothermic value for an 
asynchronous mechanism is AE — —0.375 eV, the coupling constant for synchronous transfer is J = 0.016 eV, the coupling 
constant for asynchronous transfer of weakly coupled atoms is J = 0.23 eV, the coupling constant for asynchronous transfer of 
strongly coupled atoms is J = 0.70 eV. 

above mentioned reaction of the coupled protons transfer in porphyrin molecule. Fig.|^ 24] represents the experimental 
and the theoretically estimated temperature dependencies of the reaction rate for the transfer of internal hydrogen 
(or deuteron) atoms in porphyrin. 

As it was mentioned by W. Siebrand (Div. of Chemistry, Nat. Research Council of Canada, Ottawa), a modified 
and extended microscopic model could be used to describe transfer charge transfer states in molecular dimers or 
crystals. Such states can be represented in the form A'^A~ A~A~^ that formally corresponds to a two-electron 
transfer. Two types of dimers corresponding to + and — combinations are separated by the band which is associated 
to the frequency of the transfer. None of these states have a constant dipole moment but an external field can mix 
them and induce a quite considerable dipole moment. 

In experiments, the dependence of the absorption spectrum of such systems on the external field was measured. 
The results were interpreted in terms of charge-transfer excitons. This interpretation is valid only if the external field 
is sufficiently strong, namely, to be able to break the symmetry. 

Therefore, it is of much interest to investigate splitting in the spectrum of charge-transfer dimer. In early studies 
one neglected an electron-electron interaction that appears to be of much importance here. The splitting seems to 
be much smaller than it was estimated earlier. In terms of the model presented below, electrons are attracted by 
Coulomb force of positive molecular ions and repel each other. Since the instantaneous situation when the electrons 
are at equal distance from the two molecules is energetically not favorable, a high barrier for the transfer is formed. 
An asynchronous transfer gives an intermediate state corresponding to highly excited vibronic state so that this type 
of transfer is not favorable due to Frank-Condon factors. 

Also, this state rapidly undergoes vibronic relaxation that corresponds to a non-radiative decay of charge-transfer 
states [33. As the result, it is still important to describe the observed properties in terms of controlled processes 
of electron transfer. Although W. Siebrand expressed some doubts on the possibility of a direct use of the model 
presented below to the two-proton transfer in porphyrin-type systems, he pointed out that there is a possibility to use 
such a nonlinear two-particle model in studying experimentally observed effectively two-dimensional (and two-particle) 
electron transfer. 

The microscopic model suggested in this paper allows one to interpret the experimental data, and provides complete 
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FIG. 7; Two-dimensional adiabatic potential and the path characterized by a minimal energy for (a) synchronous and (b) 
asynchronous mechanisms of NH-tautometry. 

qualitative and good quantitative description. Also, the model allows one to determine the critical point in an 
analytically exact manner (in the quasiclassical approximation), both the temperature and the coupling coefficient, 
that is an interesting theoretical result which is analogous to that of Ref . 0. As an indirect confirmation of applicability 
of the microscopic model we mention a comparison with the results of modelling of the potential energy surface in 
porphyrins [soj . 

For the case of 7-azaindole dimer, the protons are initially reside closely to nitrogen atoms in the corresponding 
5-chain "rings" of the dimer. When two coupled molecules of 7-azaindole form a dimer there may occur a cooperative 
(tunnel) two-proton transfer. The protons are moved to nitrogen atoms of 6-chain "rings" of the dimer. The proton 
transfer is accompanied by changing of the structure of orbitals of electrons. A synchronous transfer of two protons 
appears to be energetically preferable while the independent tunneling of the protons becomes forbidden due to the 
changing of the structure of orbitals. 

In Sections below, we consider a synchronous tunnel transfer of two weakly coupled particles. For two uncoupled 
particles, each of the particles moves independently in the proper double-well potential. We will study the effect of 
coupling to oscillators (an environment) at independent (asynchronous) and cooperative (synchronous) tunneling of 
the particles. Also, we will investigate fine structure of the arising two-dimensional bifurcations which occur during 
the changing of regimes. 

As it was mentioned above, the study of temperature-dependent NMR spectra of porphyrins (compounds of this 
class are known to play an important role in biological processes) revealed intermolecular regrouping (see Fig. (S)) . To 
describe this regrouping two mechanisms are suggested: (i) the synchronous one (which is associated to a synchronous 
transfer of two NH hydrogens) and (ii) the asynchronous one (with consecutive transfer of one of the particles). 

Two-dimensional potential energy surfaces (PES) V(r, tp) have been calculated by using INDO and we represent 
them in the following analytical way: 



V{r,ip) = C[a(if)(r - r, 



(6) 



The synchronous case: 



h(ip) = 1.355 + 1.681 cos(4(^)-f 0.488 cos(8(^) 
-H0.289 cos(12(/3) + 0.147 cos(16</?) + 0.049 cos(20(^) 

-h0.023 cos(24(p) eV, 



(7) 
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r™„M = 1.00- 0.14 cos(4^) A, (8) 
a((^) =38.1 cos(4(^) + 50.8 eV/A^ (9) 

The asynchronous case: 

b{ip) = 0.829 - 0.183 cos{2ip) + 0.821 cos(4(p) 

+0.260 cos(8(/3) - 0.077 cos(12(^) eV, (10) 

'^™„('^) = 1-00 - 0.14cos(4<^) 1, (11) 

a{ip) = 19.05 cos(4^) + 25.4 eV/l^. (12) 

The potential along the path characterized by a minimal energy and PES V{r^ ip) are illustrated in Fig. \7\ 

By comparing of the chosen model potential (see Fig. O with that of Fig. [3 one can see that the chosen model 
potential corresponds to a realistic PES for porphyrin. More details on PES for porphyrin have been discussed in 
Ref. [13. 



IV. TWO-PARTICLE TRANSITION PROBABILITY 



We assume that two particles independently interact with a harmonic bath. Such interaction is considered in a 
bilinear approximation. The dynamics of the environment is described by the oscillator Hamiltonian (we use h = 1, 
kg ^ 1 units with the oscillators' masses equal to one), 

i 

Each of the tunnehng particles (electrons or effective charges) interacts with the oscillator bath in the following way: 

Vi!l^{quQ.i) ^qiJ2 ^p-ph(92, Q^) - 92 E QQ„ (14) 

i i 

As in Ref. |28| . we are interested in the transition probability per unit time or, more precisely, only in its exponential 
part which can be written in the Langer's form, 

In such a consideration, metastable levels can be presented as 

r = -2Im£;, (16) 

It should be emphasized that Eq. (|15|l is valid only for temperatures below the crossover temperature while for higher 
temperatures a different prefactor is required 23] . 

Equation itT^ is obtained by generalizing the expression to nonzero temperatures 0,|l,0,IEIE0ilEIS ^> ITU 



m 111 Q 115,0^117]] 



2X;e-^°'/^Imi;, _ 2TIm^e^'^>/^ _ 2rimZ 
- e-Eo./T - ReJ2e-^'/^ ~ ReZ ' ^ ' 



Here, i labels the energy levels in the metastable state, Z is the partition function of the system, and T is the 
temperature. 

To calculate F, it is convenient to present Z in the form of a path integral |illllllli,|lllll|ll|ll|ll|ll|ll 

mniiii 

^ = n / DqiDq2DQ, exp[-5{gi, (72, QJ]. (18) 
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Here, S denotes an underbarrier action for the total system. The imaginary part ImZ is due to the decay of the 
energy states in the initial well. The validity of this approximation requires dissipation to be strong enough so that 
only an incoherent decay occurs 23] . The imaginary part in the partition function with a double well potential energy 
can be also explained due to the strong dissipation to the bath. Indeed, the particles do not come back to their 
initial state. Coherent oscillations can happen only if the interaction with bosons is weak enough jSj] or the bath is in 
nonequilibrium state [sE E3| ■ 

The integral (|18|l can be performed over phonon coordinates resulting in 



13/2 

S{qi,q2}^ J dr 



291 + 2^2+^(91,92) 



13/2 

13/2 

(19) 



+ J dT'D{r-r')[q,{T)+q2{r)][q,{r') + q2{r')]) 

-13/2 

where 



n— — oo 



(3 — h/{kBT) is an inverse temperature (below we assume that h — 1 and ks = 1), Vn — 27rn//3 is the Matsubara's 
frequency, and 

A trajectory that minimizes the Euclidean action S can be found from the equations of motion. In particular, we 
embark on the antiparallcl tunneling 

P/2 

- qi + nlqi + &iq2 + jdT'K{T - r')[9i(T') + (72 (r')] 

-13/2 

+uj^a6{-qi) - uj%e{qi) = 0, (22) 

(3/2 

- 92 + f^o92 + &iq2 + j dT'K{T ~ r')[gi(r') + 52(^0] 

-[3/2 

-L0^ae{q2) + u%e{-q2) = 0. (23) 



In Eqs. idll) and the kernel K is defined by 



-. 00 

n— — oo 

Here, is the Green's function defined by Eq. (|21|1 without the zero-frequency term, 

^(^n) = -E^+^- (25) 

Thus, we seek for sohitions of Eqs. (|22|) and (|23|1 by expanding the trajectories qi{t) and q2{t) into Fourier series, 

00 oc 

'^i - ^ E 9i^^e--, 92 = i E 9i^^e--. (26) 

n— — 00 71— — OQ 
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Introducing the renormalized frequency and interaction constant, 



fj2 (j2 

i ^ i ^ 

respectively, and substituting Eqs. H26|l into Eqs. H22() and (|23|l . we obtain that for n — 0, 

(1) , (2) 2uj^{a + b)e 
iLq + ai 

(1) (2) 2uj'^a/3 Auj'^(a + b)To 
% - 9o = 



(27) 



^0 - "1 



(28) 



and for n ^ 0, 



(1) , (2) _ 2a;^(a + b)(sin i/nTi - sin ;^„T2) 

(J) _ 2tJ^(a + 6)(sinz^„Ti + sini^„T2) 

j/„(zy4 + r2g-ai) 



Here, we have introduced the following notation: 

£ = Ti-r2, To = (ri + T2)/2. 



(29) 



(30) 



The time instants ti and T2, at which the particles pass the top points of the barrier, are determined from the following 
equations: 



gi(n) = 0, 92(7-2) = 0. 



(31) 



Equations (|H1() allow us to change the argument of the 6'-function. Namely, instead of the qi- and q2-dependencies, 
we can use a time-dependent 0-function. This reduces Eqs. H22|) and (|23|l to a linear form. 

Finally, substituting the trajectory determined from Eqs. ()26f) . (|28|l and (|29|l into Eq. (|19|) . we arrive at the following 
expression for the instanton action: 



- ai 
SmHa + bf 



E 



Sin^ VnTQ ■ COS^(t^„£/2) 



cos^ VnTQ ■ siv?{vne/2) 



{vl + ni + a^ + 2ir.)vl 



(32) 



V. PARALLEL PARTICLE TUNNELING 



For the case of parallel particle tunneling, the Euclidean action 5* can be determined similarly to antiparallel 
tunneling [see Eqs. (|22|) and H23|) ]. The trajectory minimizing the Euclidean action (instanton), can be determined 
from the equations of motion. As in the previous section, we seek for solutions of these equations in the form of the 
Fourier expansion H26I) . The time instants ri and T2 are determined by Eqs. (|31|l . 

In the case of parallel tunneling particles [the potential energy Q], the resulting Euclidean action is given as follows: 



S = 2a{a + 5)(ri + tz)^^ - -^^(a + bfin + 

u\a + bf{n-T2f 



2uj^{a + bY g 



/3 



n=l 



(w2 _ 2a)l3 
(sint^nTi + sint/„T2)^ 

(sin 1^,1 Ti - sinzy„r2)^ 



(33) 
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where ^„ is defined by Eq. if^ . 

Below, we use the fohowing notation: 



e = £*uj = (n — T2)uj, 

T = 2t*(jJ = (ti + T2)UJ, 

a* = 2a/w^ 



and assume that 



6* = 6/a, 



h> a. 



In the absence of interaction with an osciUator bath, i.e., at ^„ — 0, the action (|33|l as a function of the parameters 
e and r yields 



5 



(a + hYuj ( iar 



a + b a + b 



1-a* 



ir-\s\)a* 
1 - a* 



-coth/3* — sinh ""^/J* [cosh(/3*— r)cosh£ 
+cosh(/3*-T) - cosh(^* - |£|)] 
-(1 - a*)~3/2 [-coth(/3Vl - a*) 



+sinh" ^ (/3Vl - a*) [cosh((/3*-r)\/l - a*) 
x(cosh(£Vl - a*) - 1) +cosh((/3* - |e|)Vl - a*)]] 



As soon as the trajectory is found, Eqs. H31|l can be presented in the following form: 

sinhe[coshr coth/3* — sinhr — coth/3*] + 



1-a* 



3inh(e\/l — a*) [cosh(r\/l — a*)coth(/5*-\/l — a*) 



-sinh(TVl - a*) + coth(/3*\/l - a*)] = 
4 1 



1 + 6* 1 - a* 



1 



coshe[sinhTCoth/5* — coshr 
— 1] + sinhTCOth/3* — coshr 



+ - cosh(e-\/ 1 — a* ) [ sinh(T\/ 1 — a* ) 



xcoth(/3*Vl-a*) - cosh(TVl - a*) + 1] 
' 3inh(TVl - a*) coth(/3* %/l - a*) 



1 — a* 



-cosh(rVl - a*)] = 0. 
Simple analytic solutions of Eqs. 135(1 are obtained in the particular case when 

e = (ri - T2)lj = 0, V/3, a < 



Ti — T2 ^ — = — arcosh 
2cj 2w 



1-6* . , /3cj 

smh 

1 + 6* 2 



4' 



(34) 



(35) 



(36) 



However, a complete analysis requires numerical studies. 
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At sufficiently low temperatures, ujfi ^ 1, for 1 < 6/a < 3, and 

6 - a ^ 2a , 2(6- 
< — < a* - 



2(6 + a) - ^2 " 36 
we finally obtain, with the exponential accuracy, 

4 1 



rVT^ ~ [3 ^1(1 - a*)l/(l-vT^) 

^ 1 + 6* 1-a*^^ ^ 

1 . 4 



1 + (1 - + (3 

1 



I- a* ' 1 + 6* 

-1 



1 — a' 



-)/(i-yr^)] 



1 + 0* 1 — a* 1 — a* 



The solution (|S7|l is valid at 



/3>/3e-^^^. (38) 



We point out that an approximate solution can be found for large values of the parameter 6* (and small a*). 
However, below we focus on the more important solution H37() . The analysis indicates that there are no perturbative 
solutions of Eqs. H35|l at low temperatures and small e. 

For e = [see Eq. ^^], the action results in 



Se=o = w(6^ — a^)arcosh 
1 



1 • K^/^ 

■ smn 



2 



6* + 1 2 



cosh — 
2 



sinh-i^. (39) 



The action (|39|l coincides (up to a factor of two) with that of calculated in Ref. [23|- Thus, we have calculated the 
two-particle Euclidean action for the case of synchronous parallel motion of the two interacting particles. 
In the symmetric case (6* = 1), the action H39(l reduces to 



Se^o =4w^tanh^^. (40) 

a = b 4 

At 6* > 1, the character of the temperature dependence is almost the same. 

At e ^ 0, the corresponding action S can be obtained by substituting Eq. H37|) into Eq. (|34f) (for brevity, this 
cumbersome expression is omitted). A simple analysis shows that S^^q < 5'e=o- Moreover, it appears that the 
difference AS = Se=io ~ '5'e=o has a maximum at (jj/3 ^ 1. 

The tunnehng paths are presented by the solutions and (|S7|) . At the critical point, (3 ~ (3c, defined by 
Eq. (|38|l . a relatively abrupt change in the dynamics results in a splitting of the single trajectory (e = 0) into two 
{e ^ 0), as shown in Fig.|S| 

At (3 > (3c, i.e., for temperatures below the critical temperature, T < Tc, only a split trajectory (e ^ 0) occurs since 
Se^o<Se=o (see Fig.EJ. 

When (3 < (3c, i.e., at T > Tc, and a > ac [see Eq. (|37|) ]. the solution is determined by a single trajectory. We have 
also found that the single trajectory solution holds in the whole temperature range for a symmetric potential energy 

The above consideration is based on the analytic solutions H36f) for ti^2 near the critical point. A complete analysis 
of Eq. requires more detail. 

The numerical study of the set of transcendental equations (|35|l reveals remarkable features of two-dimensional 
tunneling. The value of the critical parameter ac decreases with temperature as shown in Figs. ^1 ^2 and 1121 The 
temperature dependence of ac indicates to the existence of a finite low-temperature limit. It becomes clear that 
the weaker the interaction between the tunneling particle, the lower temperatures are required for the synchronous 
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FIG. 8: Trajectories (a single path is characterized by e = and a splitted path is characterized by e 7^ 0) for two parallel 
moving particles, at uofi 2> 1; (l)-(4) label the projections of the minima of potential energy JHJl. 



3 



2 . 5 



1 . 5 



r 




\ 


1 


\ 




\ 


















2 









0.05 0.1 0.15 0.2 

alpha 



FIG. 9: Instanton action as a function of the interaction parameter a* — 2a /uj'^ of the two parallel moving particles, at 
ll>P S> 1; S'p = S/{uia^)\ (1) is the single trajectory; (2) is the split trajectory. 

tunneling. Since the synchronous tunneling is valid only for T < Tc [see Eq. (|38f) ] the dependence of (3c on the 
interaction parameter reveals that the minimum occurs at a* — 0.2. The nonmonotonic behavior of (3c as a function 
of 0* can be explained by the Q;*-dependence of t as well. 

Additionally, the curve in Fig. E| exhibits an anomalous behavior (an increasing part of the curve) that can be 
explained in the following way: a two-dimensional potential energy for parallel tunneling is evidently deformed by 
an increase of the interaction constant. Indeed, the minima of the potential surface become lower and the distance 
between them is larger. Small deformations in the potential energy surface change the contribution produced by the 
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FIG. 10: Numerical solutions of transcendental equations 1351 1. In addition to the studied solution t\ = T2, with larger /3 there 
are additional solutions, t\ ^ T2, shown in panels (b) and (c). Panel (b) reveals a single additional solution corresponding to 
the lower value of the Euclidean action. Panel (c) shows two different additional solutions with the one corresponding to the 
lower value of the action. 



temperature increase towards a synchronous character of the tunneling. However, for sufficiently large deformations 
of the potential, the situation is a quite opposite. The significant deformations are in favor of synchronous transfer. 
Thus, an increase in the interaction constant provides a similar effect when temperature increases. Such a behavior 
takes place up to the certain value of the interaction constant, beyond which the potential energy surface becomes 
strongly perturbed. Thus, Figs. 1111 and 1121 are complimentary to each other. Consequently, Fig. E| can be viewed as 
a bifurcation diagram. Indeed, the region below the curve corresponds to the synchronous tunneling, while the region 
above the curve corresponds to the asynchronous one. 

Fig. ^1 demonstrates how the critical parameter (3c changes from synchronous to asynchronous values with tem- 
perature. Figs. IlUb and llOb reveal the existence of the additional points, referred to as the bifurcation points, 
corresponding to the change of tunneling regimes. For the two bifurcation points shown in Fig. llUb . only one is 
physical because it corresponds to the lower value of the action (the second point is metastable). However, for suf- 
ficiently small differences between the values of the action at these bifurcation points, both close points contribute, 
thus leading to corresponding "fluctuations" in the system during the change of the regimes. For lower temperatures, 
such fluctuations become negligible because two bifurcation points contributes differently. Thus, a stable character of 
the asynchronous tunnel transfer is achieved due to a much lower value of the action at one of the bifurcation points. 

For chemical reactions mentioned in Introduction, the effect of the change in the tunneling regimes reveals as a 
cleavage in the experimental temperature dependence of the tunneling constant. Fine details of this instability have 
not been studied yet. Nevertheless, the very existence of instability at the edge of the bifurcation can be explained 
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FIG. 11; Parameter Oc as a function of inverse temperature, at the fixed value of the frequency uo and asymmetry parameter 
h*. 




FIG. 12: Parameter /3c is presented as a function of the interaction parameter of the two tunneling particles. 

as the result of specific fluctuations. A numerical analysis of the bifurcations in antiparallcl tunneling is given in the 
next section. 

VI. TWO-DIMENSIONAL ANTIPARALLEL TUNNELING 



For antiparallcl tunneling of two particles [see the potential energy Q], the instanton action as a function of the 
parameters e and r is determined by Eq. I|32|) . For ^„ — 0, we obtain 



+ 



1 - 

sinh(|e|-\/l — <5*) 
(1 -a*)3/2 



1 - (5* 



sinhjej + 



sh(eVl - a*) + 1 
(1 -a*)3/2 



x[sinh(/3Vl - a*)]-i[cosh((/3* - t)V1 - «*) 

-cosh(/?Vl - (5*)] 
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coshe — 1 
sinh/3* 



cosh(/3* - r) + cosh/3* 



(41) 



The parameters e and t are found from the following set of equations [sec Eq. H31|l ] 

— sinh£[coth/3* + coshrcoth/?* — sinhr] 



1 — a 



-sinh(e\/l - a*) [coth(/3Vl - a*) 



cosh(r\/l - (5*)coth(/3Vl ~ a*) 
+sinh(TVl - a*)] = 0, 

-1 



1 



(1 + 6*)(i-<5*) 

+ (coshe — l)(sinhrcoth/3* — coshr) 

+coshe H I [cosh(e-\/l — a*) + 11 

X [sinh(T\/l - a*)coth(/3Vl - a*) 
— cosh(T\/l — 5*)] ~cosh(e\/l ~ <5*)| = 0. 

Simple analytic solutions of Eqs. H42|l can be obtained in the following form: 

£ = (ti - T2)w = 0, V/3, 5 < 
n = r2 



(42) 



1 



arcosh 



1-6* 



— sinh( — \/\ — a* 

1 + 6* ^ 2 



T 

2w 
4' 



(43) 



2cjVl-a* 

Similarly to the case of parallel tunneling, we obtain that at low temperatures, w/? 3> 1, with the exponential 
accuracy, 

A(l -a*)i/7 



1 - (1 - a*)i/7(A/7 - (1 - <5*)-i) ^ 

1 



(44) 



Here, 



-1 



= 1 - vT^ 



(l + 6*)(l-a*) l-(5*' 

while a*, 6*, e, and r are determined as for the parallel transfer. 

The solution (|44|l is valid at a*^^ < a* < (5*2 , where the lower and upper bounds a*]^ and 5*2 are derived from 
a cumbersome transcendental equation (for brevity it is not presented here). Particularly, in the symmetric case, 
6* = 1, we obtain the condition in a simple analytic form, 1/4 < 2a /uP' < 1. 

Furthermore, an approximate solution can be found for large values of the parameter 6* — b/a (and small a*). 
However, we restrict our analysis to the more important physical solution H44(l . 

The /3-dependent solution (|44|l is valid for (3 > (3c, where 



In 



U/7 



l_(l_«*)l/7(^_^)- 



(45) 



At uj(3 ^ 1, the solutions of Eqs. I|42|) can be found perturbatively (for small e) with given values of the parameters 
(6 - a)/(6 + a) and a*. At e = [the solution (03))], the action iglj yields: 



S 



uip - a2) 
(1 -5*)3/2 



arcosh 



a . LoP^/l — a* 
— smh 



bj{h + a) 



2 r 



2(1-52) (l_5*)3/2 



coth 



rsinh-^^^MZ^ 



(6 + ay 



(46) 
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FIG. 13: Trajectories (the basic path characterized by e = and the spHt one is characterized by e 7^ 0) at Luf3 S> 1 of two 
antiparallel tunneling particles; (l)-(4) denote the projections of the minima of the potential energy surface Ua,{qi,q2) defined 
by Eq. |0J . 



For the symmetric potential (a ~ b) and £ = 0, we obtain that (see Fig. Ill|l 

^=(1— 4 ■ 

We do not present here a cumbersome expression for S^^o which one can obtain by substituting the solutions r 
and e into Eq. H41I) . A simple analysis reveals that S^^q > Se=o for /3 > (3c and for relevant a* . Similarly to parallel 
transfer, the tunnel paths can be found from Eqs. H43|) and 144|) . These trajectories on the {qi, q2)-pl&iie are shown in 

Fig. ini 

As for parallel tunneling and at /3 > /3c, the pair-tunneling changes from a single to double trajectory regime. In 
contrast to the parallel tunneling, such a splitting occurs at any values of the parameters of the potential. At (3 > (3c, 
we have S^^q > S^^o and, therefore, 5^=0 determines the tunneling rate. At /3 < /3c, the two degenerated trajectories 
are transformed into a single trajectory, qi = —(72, corresponding to synchronous antiparallel transfer. 

For single particle tunneling, there is only a single tunneling path (instanton) minimizing the action. Hence, there 
are two different types of trajectories for the pair of interacting particles. Namely, the main contribution to the 
instanton action is determined either by the single or by the double-degenerated paths depending on the value of /3. 
We also point out that in the case of the parallel tunneling, the particles do not simultaneously pass the top points 
of the barrier, ri =/= T2. This means that the tunneling transfer is asynchronous. 

At small values of the interaction parameter a* [see Eq. and at temperatures such that (3 < (3c [see Eq. 
there is no splitting of a single path {qi — (72)- Therefore, the particles pass the top of the barriers at the same 
instants (ri = T2). Consequently, the transfer of the particles is synchronous. The temperature dependence for the 
antiparallel transfer action is plotted in Fig. ^] at various a*. 

The type of the interaction given by Eqs. 101-© is such that it does not affect the motion along the "center of 
mass" coordinate, qi = q2. For this reason, the Euclidean action is independent of the interaction parameter as for 
parallel transfer. Since the state of the interacting system characterized by a maximal value of the relative coordinate, 
qi = ~'q2, is preferable (as it provides the lower action), it becomes clear that the instanton action decreases with 



19 



2.5 

2 

^ 1.5 



1 

0.5 



1 2 3 4 5 6 

w 



FIG. 14: Instanton action for antiparallel tunneling (e = 0, a = 6) as a function of an inverse temperature at the two different 
values of the interaction parameter; Sa ~ S/i/^Loa?), w = lij/3/4; (1) a* — 0.5; (2) 3* = 0.1; (3) a dashed hue corresponds to 
the action 14011 for the parallel transition (e = 0, a = b). 

the interaction parameter in the paraUel transfer along the degenerated tunnel trajectories, and increases with the 
interaction parameter for the antiparallel tunneling. 

For antiparallel tunneling, synchronous transfer (ti — T2) takes place, while asynchronous transfer is forbidden due 
to the greater contribution to the Euclidean action (see Fig. 1141) . 

The validity condition for weakly interacting instanton-antiinstanton pairs, in the adiabatic approximation, was 
discussed in Ref. j27]. 

The numerical analysis of transcendental equation H42|l reveals interesting features for a transition region between 
the tunneling regimes, i.e., a fine structure near the first bifurcation point for the antiparallel transfer. The numerical 
results are presented in Fig. ^| We found that, in addition to the first bifurcation point characterized by the two 
solutions (Fig. Ell), there exist additional bifurcation points at lower temperatures, i.e., the four pairs fFig.lTSb'). the 
six pairs fFig.llSb). and there are even twelve pairs of additional solutions at /3* = 19.2009 {a* — 0.05), etc. We refer 
to this phenomenon as multiplication of bifurcations or a cascade of bifurcations. Such an effect resembles a scenario 
of transition to chaos. 

Although the synchronous regime is preferred due to the minimal instanton action, in a certain temperature range 
its value is comparable with those of corresponding to the cascade solutions. As a result, quantum fluctuations of a 
nonregular character occur in contrast to the parallel transfer. The antiparallel tunneling is, thus, characterized by 
the instability of the transition due to the synchronous to asynchronous behavior. Such instabilities are similar to 
continuous second-order phase transition, while the parallel tunneling is viewed as a step process similar to phase 
transition of a first order (see Fig.llOfl. The dependencies shown in Figs.lllland ll2l /3c(a) and ctdP), for the antiparallel 
transfer are found to be of the same character as those of the parallel transfer. 

In summary, our investigation reveals a quite complicated fine structure of the transition for parallel and antiparallel 
tunneling of two particles with different degenerate trajectories leading to the bifurcation cascade. 

VII. EFFECT OF A PROMOTING MODE 

In this section, we study the effect of a heat bath on tunneling transition of two interacting particles. In many 
tunneling reactions, interaction with a vibrational subsystem can often be approximated by interaction with a single 
vibrational mode (a so-called promoting mode). As it follows from Eq. H19|l . a heat bath affects only the dynamics 
of the center of mass {qi =92)- Therefore, in the case of an antiparallel motion, the medium does not affect the rate 
constant, while for parallel tunneling it essentially contributes to the transfer rate. For both parallel and antiparallel 
transfer, bilinear interaction of the particles with a single oscillator can make a qualitative change in the character of 
tunneling. 

At small values of the interaction parameter between the two tunneling particles dissipative effects become important 
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FIG. 15: Numerical solutions of transcendental equation 14211 . In addition to the studied solution ri — T2, with respect to 
/3, there are the additional solutions, ri T2, shown in panels (a), (b) and (c), that correspond to one, four, and six (pairs) 
additional solutions, respectively. 
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|27| . In two dimensions the dissipative effects are more pronounced for parallel rather than antiparallel tunneling. 
The latter increases with temperature with a considerable contribution to the prefactor. In the present work, we 
are interested in the tunneling rate assuming only exponential evolution of the transition probability. However, 
nonexponential evolution can occur in an noncquilibrium environment [sM l39l l40l . Such a case is not discussed 
here. Accordingly, a reservoir is assumed to be in thermodynamic equilibrium, i.e., the tunneling transition is rather 
slow compared to the thermodynamic relaxation. Thus, we assume that dissipation affects the value of instanton 
action only. 

For the case of the antiparallel tunnel transfer, the action H32|l can be calculated with the following vibronic Green's 
function: 



D{v^) = — 



C2 



(48) 



where Wi, is the frequency of the vibrational mode. After some tedious calculations, one obtains the following 
expression for the instanton action H32|): 



S = 



"0 



/3sinh(r>o(/3/2-ro))sinh(r>oro) 
2173sinh(17o/3/2) 



A /3(02 - w2 )cosh(r2,(/3/2 - To))cosh(r2,;To) 
/3e (3e /3sinh(£ilo) 



4(t^2_2C2/^i) 



ml 



ml 



P{nl - LuDsinhien,) /?cosh(er2o) 



^ {-lyAn^inj-nl) 4173sinh(/317o/2) 
x[cosh(r2o(/3/2 - 2to)) - cosh(r2o/3/2)] 



(49) 



^ {-lymiini - r22)sinh(/3rj,/2) 

X [cosh(a,(/3/2 - 2ro)) + cosh(f7i/3/2)] } . 
Here, we have introduced the following notation: 

nl^uj^- a^ nl = i(c.2+^i + ^(c.2-^2)2+8c^2 )^ 

For the parallel tunneling transition, a corresponding action can be found in a similar way. The action H34I) as a 
function of the parameters e* = ri — T2 and r* = (ti + T2)/2, with the vibronic frequency ojl and coupling constant 
C, yields (see also Ref. "27^): 



S={a+b){3a-b)uj\*- 



2^*2 



2(cj2_2q) 13 
Aiu^ia + bfs*^ uj^a + bf ^ [-Ifiuj^ - ^3-^) 



(tj2 _ 2a)(3 
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coth(^) - sinh-i(^) (cosh[(^ - 2r*)^ 



cosh[(^ - e*)^/x^ 



-cosh[( — 
2 ^4 



^icosh[(^-2r*+e*)V^ 



-coth( ^ ^uj-^-2a) + sinh" ^ ( ^ V - 2a) 



2T*-e*)V^ 

uo^{a+bf 
' 2(t^2_2a)3/2 
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^-cosh[(^-2r*)v/^;^2i:2^] + cosh[(^-e*)v/^j2Z^] 

+ icosh[(^ - 2r* - e*)V^2_2a] 
+ icosli[(| - 2t* + £*)v/c^2_2a]^ 



(50) 



Here, we have denoted 



1. 2 2 1~ 

/ C2 
7= J(^2+^2 +_)2_4^2^2^ 

V 

For particular values of the interaction constant a and in the absence of interaction with the oscillator bath, the critical 
temperature Tc (at which the synchronous and asynchronous tunnel regimes interchange) is found from Eqs. (|38|l and 
(|45|l . These equations can be generalized to non-zero interaction with the promoting mode. Typically, the critical 
temperature is found to be in the range from 10 K to 400 K. In glasses, Tc can be very small while for chemical 
reactions it can be rather large. Additionally, Tc depends on the mean distance between the particles and, therefore, 
on their concentration. 

Quantum tunneling is important ^28] when ksTc/ (huj) < 1. Therefore, the symmetry breaking effects can take 
place at relatively high temperatures depending on the "frequency" of a barrier. For example, for porphyrins the 
critical temperature Tc is estimated to be about 200 K. 



VIII. CONCLUSIONS 



In the single instanton approximation, we have calculated the Euclidean action (|19|l for the models characterized 
by the different adiabatic potential energy surfaces, Q and (jsj, and made a detailed comparative analysis of the 
tunneling rate for two interacting particles moving in parallel or anti par allel within a dissipative environment. We 
have also assumed exponential evolution of the transition probability 123, |33, l39l l40l l4ll | . 

We have shown that the change in a tunneling regime from synchronous to asynchronous transfer for parallel 
transition occurs as a step process, similar to phase transition of a first order, while for the antiparallel transfer it 
resembles second order phase transition. 

We have explained the effect of a cleavage in the experimentally observed [23. l25l |2^ l33l| temperature dependence 
of the reaction rate for two tunneling particles. It has been shown that the effect of symmetry breaking is stable for 
the parallel and unstable for the antiparallel transfer, as it is observed experimentally for some porphyrin systems 
[23. 125I I29. l3a | . We have found a complicated fine structure in the bifurcation region due to the quantum fluctuations 
for the parallel two-dimensional tunneling. For the antiparallel tunneling, the contribution of four, six, twelve, etc., 
pairs of trajectory becomes important resembling the transition to chaos. 

Additionally, we have studied interaction of two particles with phonons. Such coupling essentially modifies 
the antiparallel and parallel transitions in different ways. As it follows from Eq. (|19|l . the interaction with the 
reservoir does not change the dynamics of the center of mass for the antiparallel motion, while it makes a significant 
contribution to the transfer rate for the parallel transfer. Finally, Eq. (|38|l determines the validity condition for 
temperatures beyond of which stable two-dimensional synchronous tunneling correlations of all the kinds occur. 

The authors would like to thank A. I. Larkin and B. I. Ivlev for the stimulating interest in this work. 
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